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We study the Casimir problem as the limit of a conventional quantum field theory coupled to a 
smooth background. The Casimir energy diverges in the limit that the background forces the field 
to vanish on a surface. We show that this divergence cannot be absorbed into a renormalization of 
the parameters of the theory. As a result, the Casimir energy of the surface and other quantities like 
the surface tension, which are obtained by deforming the surface, cannot be defined independently 
of the details of the coupling between the field and the matter on the surface. In contrast, the energy 
density away from the surface and the force between rigid surfaces are finite and independent of 
these complications. 
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The vacuum energy of fluctuating quantum fields that 
are subject to boundary conditions has been studied in- 
tensely over the half-century since Casimir predicted a 
force between grounded metal platesP,|3,l3- The plates 
change the zero-point energies of fluctuating fields and 
thereby give rise to forces between the rigid bodies or 
stresses on isolated surfaces. The Casimir force between 
grounded metal plates has now been measured quite ac- 
curately and agrees with his prediction 

Casimir forces arise from interactions between the fluc- 
tuating fields and matter. Nevertheless, it is traditional 
to study idealized "Casimir problems" where the physi- 
cal interactions are replaced ah initio by boundary condi- 
tions. In this Letter we study under what circumstances 
this replacement is justified. A real material cannot con- 
strain modes of the field with wavelengths much smaller 
than the typical length scale of its interactions. In con- 
trast, a boundary condition constrains all modes. The 
sum over zero point energies is highly divergent in the 
ultraviolet and these divergences depend on the bound- 
ary conditions. Subtraction of the vacuum energy in the 
absence of boundaries only removes the worst divergence 
(quartic in three space dimensions). 

The fact that the energy of a fluctuating field diverges 
when a boundary condition is imposed has been known 
for many years|3, Schemes have been proposed to 
cancel these divergences by introducing new, ad hoc sur- 
face dependent counterterms '9J or quantum boundary 
functions [injl. We are not interested in such a formal 
solution to the problem. The method of renormalization 



in continuum quantum field theory without boundaries 
(QFT) provides the only physical way to regulate, dis- 
cuss, and eventually remove divergences. Therefore we 
propose to embed the Casimir calculation in QFT and 
study its renormalization. After renormalization, if any 
quantity is still infinite in the presence of the boundary 
condition, it will depend in detail on the properties of 
the material that provides the physical ultraviolet cutoff 
and will not exist in the idealized Casimir problem. Sim- 
ilar subtleties were addressed in the context of dispersive 
media in Ref. 0| . 

It is straightforward to write down a QFT describing 
the interaction of the fluctuating field (p with a static 
background field cr(x) and to choose a limit involving 
the shape of tT(x) and the coupling strength between 
<j) and CT that produces the desired boundary conditions 
on specified surfaces. We have developed the formal- 
ism required to compute the resulting vacuum energy in 
Ref. UJ. Here we focus on Dirichlet boundary condi- 
tions on a scalar field. Our methods can be generalized 
to the physically interesting case of conducting boundary 
conditions on a gauge field. 

Ideally, we seek a Casimir energy that refiects only the 
effects of the boundary conditions and not on any other 
features of cr(x). Therefore we do not specify any action 
for a except for the standard counterterms induced by 
the <j)-a interaction. The coefficients of the counterterms 
are fixed by renormalization conditions applied to per- 
turbative Green's functions. Once the renormalization 
conditions have been fixed by the definitions of the phys- 
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ical parameters of the theory, there is no ambiguity and 
no further freedom to make subtractions. Moreover, the 
renormahzation conditions are independent of the par- 
ticular choice of background (t(x), so it makes sense to 
compare results for different choices of cr(x), ie. differ- 
ent geometries. Having been fixed in perturbation the- 
ory, the counterterms are fixed once and for all and must 
serve to remove the divergences that arise for any phys- 
ically sensible (t(x). The Yukawa theory with coupling 
g in three space dimensions gives a textbook example: 
The gtpcFip coupling generates divergences in low order 
Feynman diagrams proportional to a^, cr'* and {da)^ and 
therefore requires one to introduce a mass, a quartic self- 
coupling, and a kinetic term for a. This is the only con- 
text in which one can study the fluctuations of a fermion 
coupled to a scalar background in three dimensions. 

In this Letter we study the vacuum fluctuations of a 
real scalar field (j) coupled to a scalar background cr(x) 
with coupling A, £int(</>, cr) = iAcr(x) 0^(x, i). In the 
limit where (t(x) becomes a (5-function on some surface S 
and where A ^ cx), it is easy to verify that all modes of 
(j) must vanish on S. We call this the Dirichlet limit. It 
consists of the sharp limit, where cr(x) gets concentrated 
on S, followed by the strong coupling limit, A — > cx). In 
general, we find that the divergence of the vacuum energy 
in the Dirichlet limit cannot be renormalized. Generally, 
even the sharp limit docs not lead to a finite Casimir 
energy except in one dimension, where the sharp limit 
exists but the Casimir energy diverges as A In A in the 
strong coupling limit. 

This divergence indicates that the Casimir energy of a 
scalar field forced to vanish on a surface in any dimension 
is infinite. However, all is not lost. The unrenormalizable 
divergences are localized on S, so quantities that do not 
probe S are well defined. For example, it is straightfor- 
ward to show that the vacuum energy density away from 
S is well defined in the Dirichlet limit, even though the 
energy density on S diverges Jl^. We expect that this is 
true in general. The forces between rigid bodies are also 
finite in the Dirichlet limit. But any quantity whose def- 
inition requires a deformation or change in area of S will 
pick up an infinite contribution from the surface energy 
density and therefore diverge. For example, we will see 
explicitly that the vacuum contribution to the stress on a 
(generalized) Dirichlet sphere in two or more dimensions 
is infinite. 

The remainder of this Letter is organized as follows: 
First we briefly review our computational method and 
discuss the structure of the counterterms required by 



renormahzation. Then we present two examples, leav- 
ing the details to Ref. We begin with the simplest 
Casimir problem: two Dirichlet points on a line, where 
we can compare our results with standard calculations 
that assume boundary conditions from the start 2] . We 
find that the renormalized Casimir energy is infinite but 
the Casimir force is finite in the Dirichlet limit. We show 
how the QFT approach resolves inconsistencies in the 
standard calculation. Next we study the Dirichlet circle 
in two dimensions. We demonstrate explicitly that the 
renormalized Casimir stress on the circle diverges in the 
sharp limit. We show that the divergence is associated 
with a simple Feynman diagram and will persist in three 
dimensions (the "Dirichlet sphere") and beyond. 

We define the bare Casimir energy to be the vac- 
uum energy of a quantum field (j) coupled to a back- 
ground field a by £i,it(0, cr), minus the vacuum energy 
in the absence of a. This quantity can be written as 
the sum over the shift in the zero-point energies of all 
the modes of </) relative to the trivial background a — 0, 
Ehare[o-] = tZ]n('^"H -^n^)- Equivalcntly, using the 
effective action formalism of QFT, £'baio[o'] is given by 
the sum of all 1-loop Feynman diagrams with at least 
one external a field. The entire Lagrangian is 

jr ^ id.cpd^cj^ - - ^</.V(x) + /:ctM , (1) 

where Cqt [c] is the counterterm Lagrangian required by 
renormahzation. Combining its contribution to the en- 
ergy with -Ebarc [c] yields the renormalized energy Ecas [c] ■ 
We have taken the dynamics of the background field 
(t(x) to include only the 4>-a coupling and the countert- 
erms required by renormahzation theory. The frequencies 
{ci;[tT]} are determined by —V^0(x)-f(m^-fAcr(x))0(x) = 
tj2[cr]0(x). This is a renormalizable quantum field the- 
ory, so i?cas[c] will be finite for any smooth a and fi- 
nite A. We use the method developed in Ref. to 
compute the Casimir energy of the background config- 
uration exactly while still performing all the necessary 
renormalizations in the perturbative sector. The inter- 
ested reader should consult Ref. 1^ for an introduction 
to the method and Ref. 15] for applications. We assume 
that the background field cr(x) is sufficiently symmetric 
to allow the scattering amplitude to be expanded in par- 
tial waves, which we label by £. We express the renormal- 
ized Casimir energy as a sum over bound states ujij plus 
an integral over continuum modes with uj = \/W^^^^rr? , 
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where Ni denotes the multiplicity, 5i the scattering phase 
the continuum density of states in the 



shift and 



1 dSt_ 
dk 



partial wave. The subscript N on indicates that 
the first iV terms in the Born expansion of 8i have been 
subtracted. These subtractions are compensated exactly 
by the contribution of the first TV Feynman diagrams, 

■^FD = E^i ^FD- III eq. lEJ i^cT is the contribution of 
the counterterm Lagrangian, £ct- Both and i?cT 
depend on the ultraviolet cutoff 1/e, but + Eq^ re- 
mains finite as e — !■ 0. One can think of e as the stan- 
dard regulator of dimensional regularization, although 
our methods are not wedded to any particular regular- 
ization scheme. After subtraction, the fc-integration in 
eq. (121) converges and can be performed numerically for 
any choice of cr(x). It is convenient for computations to 
rotate the integration contour to the imaginary fc-axis 
giving 



dt 



t 



(3) 

where t = —ik. The real function I3g{t) is the log- 
arithm of the Jost function for imaginary momenta, 
I3i{t) = InF^(ji). Efficient methods to compute (3i{t) 
and its Born series can be found in Ref. J^j . The renor- 
malized Casimir energy density for finite A, ecas(x), can 
also be written as a Born subtracted integral along the 
imaginary fc-axis plus contributions from counterterms 
and low order Feynman diagrams [l^ . 

In less than three dimensions only the lowest order 
Feynman diagram diverges, so only a counterterm lin- 
ear in (7 is necessary, £ct = ciAcr(x). Since the tadpole 
graph is also local, we can fix the coefficient ci by re- 
quiring a complete cancellation, i?pp -I- Ect — 0. In 
three dimensions it is necessary to subtract two terms in 
the Born expansion of [3i{t) and add back the two low- 
est order Feynman graphs explicitly. The counterterm 
Lagrangian must be expanded to include a term propor- 
tional to cr^, £cT = ciAct(x)-|-C2^ cr^(x). The new term 
cancels the divergence generated by the vacuum polariza- 
tion diagram -Ep^, but it does not completely cancel -Epo, 
because E^^ it is not simply proportional to J d?xa^{^). 

To fix C2 we can only demand that it cancels E^^ at a 
specified momentum scale = A'P. Different choices of 
M correspond to different models for the self-interactions 
of cr and give rise to finite changes in the Casimir energy. 

We use eq. (jSJ and the analogous expression for the en- 
ergy density with backgrounds that are strongly localized 
about S, but not singular, to see how the Dirichlet limit 
is approached. It is straightforward to relate the Casimir 
energy density at the point x to the Green's function at 
X in the background cr, and then to show that it is fi- 
nite as long as cr(x) = 0. Thus we find that the Casimir 
energy density at any point away from S goes to a fi- 



nite limit as cr — > ^^(x) and A ^ oo and that the result 
coincides with that found in boundary condition calcula- 
tions. We also find a finite and unambiguous expression 
for the renormalized Casimir energy density where cr(x) 
is nonzero, as long as it is nonsingular and the coupling 
strength is finite. But as we approach the sharp limit, 
the renormalized energy density on S diverges, and this 
divergence cannot be renormalized. 

By analyzing the Feynman diagrams that contribute 
to the effective energy we can deduce some general re- 
sults about possible divergences in the Casimir energy 
and energy density in the sharp limit. In particular, the 
divergences that occur in the Casimir energy in the sharp 
limit come from low-order Feynman diagrams. Specif- 
ically, using dimensional analysis it is possible to show 
that in n space dimensions the Feynman diagram with 
m external insertions of cr is finite in the sharp limit if 
m > n. 

Although less sophisticated methods can be used to 
obtain the energy density at points away from S where 
renormalization is unnecessary, as far as we know only 
our method can be used to define and study the Casimir 
energy density where cr(x) is nonzero and therefore on S 
in the Dirichlet limit. 

Consider, as a pedagogical example, a real, massive 
scalar field (l)(t, x) in one dimension, constrained to vanish 
at X ~ —a and a. The standard approach, in which 
the boundary conditions arc imposed a priori, gives an 
energy 0] 
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where the tilde denotes the imposition of the Dirichlet 
boundary condition at the outset. From this expression 
one obtains an attractive force between the two Dirichlet 
points, given by 



F(a) 



dE2 
' d{2a) 



dt 
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TT ^t^ - m2(e4"* - 1) 



(5) 



In the massless limit, we have E2{a) = — 7r/48a and 
F{a) = -7r/96a2. 

These results are not internally consistent, suggesting 
that the calculation has been oversimplified: As a — > oo, 
£■2 (a) — * —TTi/'^^ indicating that the energy of an isolated 
"Dirichlet point" is —to/4. As a —+ we also have a 
single Dirichlet point, but E2{a) — > oo as a — > 0. Also 
note that i?2(a) is well defined as m —^ 0, but we know on 
general grounds that scalar field theory becomes infrared 
divergent in one dimension when m —^ 0. 

We study this problem by coupling <j){t, x) to the static 
background field cr(x) = 6{x + a) + d{x — a) with coupling 
strength A as in eq. . An elementrary calculation gives 
the renormalized Casimir energy for finite A, 



The same method can be applied to an isolated point 



givmg, 
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For any finite coupling A, the inconsistencies noted in 
E2{a) do not afflict i?2(a, A): As a ^ oo, E2{a,X) 
2Ei(\), and as a ^ 0, E2(a,\) -> Ei{2\). Also £^2(0, A) 
diverges logarithmically in the limit m as it should. 
The force, obtained by differentiating eq. |(BJ with respect 
to 2a, agrees with eq. (jSJ in the limit A — > 00. 



Note, however, that £'2(0, A) diverges like A log A as 
A — !■ 00. Thus the renormalized Casimir energy in a sharp 
background diverges as the Dirichlet boundary condition 
is imposed, a physical effect which is missed if the bound- 
ary condition is applied at the outset. 



The Casimir energy density for x ^ ±a can be calcu- 
lated assuming Dirichlet boundary conditions from the 
start simply by subtracting the density in the absence of 
boundaries without encountering any further divergences 
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e2(x, a) = 
e2(x, a) = 



dtVt^ 



9 00 r ri.TT 

\ - cos ' 



HL(x-a)] 



4a 



777- 

KQ{2m\x — a\) for \x\ > a . 

27r 



for \x\ < a 



(8) 



The Casimir energy density for finite A was computed 
in Ref. In the limit A ^ 00 it agrees with eq. © 

except at X = ±a where it contains an extra, singular 
contribution. If one integrates eq. over all x, ignoring 
the singularities at a; = ±a, one obtains eq. Q). Including 
the contributions at ±a gives eq. 

This simple example illustrates our principal results: 
In the Dirichlet limit the renormalized Casimir energy 
diverges because the energy density on the "surface," x — 
±a diverges. However the Casimir force and the Casimir 
energy density for all x ^ ±a remain finite and equal to 
the results obtained by imposing the boundary conditions 
a priori, eqs. (O and 

A scalar field in two dimensions constrained to vanish 
on a circle of radius a presents a more complex prob- 
lem. We decompose the energy density in a shell of 
width dr at a radius r into a sum over angular momenta, 
= X^fco ^^('')' where ee{r) can be written as an in- 
tegral over imaginary momentum t — —ik of the partial 
wave Green's function at coincident points Ge{r, r; it) and 
its radial derivatives. First suppose we fix (t(x) = S{r—a) 
and consider r ^ a. It is easy to see that the difference 
\Gi(r, r, ji)]o between the full Green's function (r, r, it) 
and the free Green's function cf^ (r, r, it) vanishes expo- 
nentially as < — *■ 00. For finite A, both the t-integral and 
the £-sum are uniformly convergent so A — > 00 can be 
taken under the sum and integral. The resulting energy 
density, e(r), agrees with that obtained when the Dirich- 



let boundary condition, 4>{a) = 0, is assumed from the 
start. As in one dimension, nothing can be said about 
the total energy because e(r) is not defined at r = a, but 
unlike that case the integral of e(r) now diverges even in 
the sharp limit for finite A. 

To understand the situation better, we take ct(x) to 
be a narrow GauBian of width w centered at r = a 
and explore the sharp limit where w — > and (t(x) — > 
(5(r — a). For u) 7^ 0, cr does not vanish at any value 
of r, so [G£(r, r, it)]o no longer falls exponentially at 
large t, and subtraction of the first Born approxima- 
tion to Gi{r,r,it) is necessary. As in one dimension, 
the compensating tadpole graph can be canceled against 
the counterterm, ciAcr(x). The result is a renormal- 
ized Casimir energy density, e(r, w. A), and Casimir en- 
ergy, E(w, A) = Jq dre{r, w, A), both of which are finite. 
However as w ^ both e(a, w, A) and E{w, A) diverge, 
indicating that the renormalized Casimir energy of the 
Dirichlet circle is infinite. 

The divergence originates in the order A^ Feynman di- 
agram. We study this diagram by subtracting the second 
Born approximation to Gi{r,r,it) and adding back the 
equivalent diagram explicitly. Then the ^-sum and t- 
integral no longer diverge in the sharp limit. In the limit 
w —> 0, the diagram contributes 

firn^ Ep^ — — — J dp Jg {ap) arctan ^ (9) 
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which diverges Hke In A. The divergence originates in the 
high momentum components in the Fourier transform of 
<j{r) — S(r — a), not in the high energy behavior of a foop 
integral, and therefore cannot be renormahzed. Taking 
X —^ oo only makes the divergence worse. Because it 
varies with the radius of the circle, this divergence gives 
an infinite contribution to the surface tension. This diver- 
gence only gets worse in higher dimensions (in contrast to 
the claim of Ref. For example, for a{r) — 6{r — a) 

in three dimensions the renormalized two point function 
is proportional to an integral over p of a function pro- 
portional to A^o^p^Jq (pa) Inp at large p. The integral 
diverges like A In A. Such divergences cancel when we 
compute the force between rigid bodies, but not in the 
case of stresses on isolated surfaces. 

In summary, by implementing a boundary condition 
as the limit of a less singular background, we are able to 
study the divergences that arise when a quantum field is 
forced to vanish on a prescribed surface. For all cases we 
have studied, the renormalized Casimir energy, defined 
in the usual sense of a continuum quantum field theory, 
diverges in the Dirichlet limit. Physical cutoffs (like the 
plasma frequency in a conductor) regulate these diver- 
gences, which are localized on the surface. On the other 
hand the energy density away from the surfaces or quan- 
tities like the force between rigid bodies, for which the 
surfaces can be held fixed, are finite and independent of 
the cutoffs. Observables that require a deformation or 
change in area of S cannot be defined independently of 
the other material stresses that characterize the system. 
Similar studies are underway for fluctuating fermion and 
gauge fields, leading to Neumann and mixed boundary 
conditions with the same types of divergences. 
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